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Abstract 

The ratio of average muhiphcities in gluon and quark jets is shown 
to become noticeably smaher in higher-order QCD compared to its 
lowest order value what improves agreement with experiment. QCD 
anomalous dimension has been calculated. It has been used to get 
energy dependence of mean multiplicities. 
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The ratio of average multiplicities in gluon ( (nc) ) and quark ( {np) ) 
jets ( let us call it r = (nc) / {np) ) was of much debate during last years. 
Its value Nc/Cp = 9/4 predicted in the lowest order approximation of QCD 
is too large to fit experimental data. The simplest correction |jl], |[ has 
reduced this value to 2.05. Recently it has been shown 0] that the exact 
solution of fixed coupling QCD equations gives rise to ever smaller value 
r = 1.84 ± 0.02. The above values refer to the ratio of the number of partons 
in jets. In Monte Carlo computation using definite hadronization schemes 
they are recalculated to hadronic ratios observed in experiment which are 
about 40 percents smaller. Using one of such models ( namely, Herwig), it 
has been shown in that the above partonic value of r agrees quite well 
with experimental hadronic value 1.27 ± 0.04 ± 0.06. 

We show here that quite similar values are obtained for the partonic ratio 
in the running coupling QCD if higher order terms are taken into account. 
Therefore, one concludes that this ratio is rather insensitive whether the 
coupling constant is fixed or running but requires proper treatment of energy 
conservation. 

We have calculated also the QCD anomalous dimension in the same ap- 
proximation and have found the energy dependence of average multiplicity. 
They are more sensitive to the running property of the coupling constant. 

Let us remind briefly main notations. The generating functions of gluon 
Gciy, z) and quark Gp{y, z) jets are 

oo oo 

Gaiy^z) = J2il + zrPniy), Gp{y,z) = Y.{1 + zT {y) , (1) 

ra=0 n=0 

where Pn{y) and Pn{y) are the parton multiplicity distributions in gluon 
and quark jets, correspondingly, y = ln(Q/Qo) ? Q is the jet virtuality, 
Qq =const. The normalized factorial moments Fg and $g are related to 
generating functions as: 

oo g °° 

GG = E-j{ncyFg, Gp = j:-{npr%, (2) 

q=0 y- q=0 

The generating functions satisfy the following integro-differential equations 
i 
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G'a{y) = C dxK^{x)^l[GG{y + \nx)GG{y + ln(l - x)) - Gdy)] 
Jo 

+nf f dxK^{x)^l[GF{y + lnx)Gf (y + ln(l - x)) - Gaiy)], (3) 

G'p{y)= I' dxK${x)^l[GG{y + \nx)GF{y + \n{l-x))-GF{y)l (4) 
Jo 

where G'{y) = dG/dy, 7q = 2Ncas/'n', o^s = 27r/(ll — \nf)y is the running 
couphng constant, n/ is the number of active flavours of quarks. 
The kernels of equations are given by QCD as : 
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[x^ + {l-xY], 
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(5) 
(6) 
(7) 



where iVc = 3 is the number of colours and Gp = {N^ — l)/2Nc = 4/3. 

In solving eqs.(^,(^) we follow the method developed in 0, |^ and expand 
the generating functions in Taylor series at large y. We leave aside all terms 
with second and higher derivatives here because their contribution to average 
multiplicities is of higher order in coupling constant as can be easily checked. 
In that case, eqs.(|D,@ become 



12N, 



{\nGG{y)r= / dy'j',{y')[GGiy')-l] 
Jo 

-/ii7o[G'g(i/)-1] + MoG"c(i/) 



Gciy) 



- 1 



GG{y) 



(In GHz/))' , 



(8) 



(lnGHl/))'=^ 



dy'i',{y')[GG{y') 



-\{GG{y) - 1) + \G'G{y) - h,{\nGF{y))'GG[y) 



(9) 
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where hi = 11/24, /i2 = (67 - 67r2)/36, h = {An^ - 15)/24, = 13/3. 

To calculate average multiplicities, one should keep only the terms with 
q = and 1 in (^. Normalization condition is Fg = Fi = $o = = 1 • 
After differentiation, one gets from 



{nc)" = 7o {(1 + ^Bha'o){nG) - 2h{nGy + hinc)" 
i^[2{2M'-{nay)-K{n,)" 

-4572(2(n^)-(nG))]} 



UN, 



(10) 



{ufY 



^lo {(1 + Is^oXnc) - line)' + line)" - h,{n,)"] , (11) 



where 7o = -B%, B = {fN^ - ^nf)/8Nc has been substituted. The QCD 
anomalous dimension ■y{y) is defined as 

ry 

{no) = exp( J -f{y')dy'). 
We represent 7 and r expanding them in 70 as 

7 = 7o(l - ai7o - a27o) + C'(7o)> 



^ = ^(1-^70 -r27o) + 0(7o')- 



(12) 

(13) 
(14) 



Since 7q ~ 79, we get r' = up to terms 0(7q). However, such terms are 
taken into account in the formulae below. Up to correction terms of order 
0(7o) one gets 



{no) 



7o 



1 - ai7o - a27o 



1 - (2ai + 5)70 - (2a2 -a\- 2aiB)% 



(15) 
(16) 



1 + (ri - ai)7o + (r2 + - a2 - - Bn)-/^ , (17) 
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l + in-B- 2ai)7o + (r2 + + a?- 
3Bri - 2airi - 2a2 + 2Bai)^^ 



(18) 



We insert (|Tll)-(0) into (|T0|),(|TT1) and keep again the terms of order 0(7q) 
Therefrom we obtain after simple algebraic manipulations 



ai = hi + 



12m 



B 
~2'' 



3 
4' 



(19) 
(20) 



a2 = --{h2 - al - 2ai \B - [hi 



12N} 



UN, 



[hi + 3 - iB] + 4:B[hi + 



n. 



12iVf 



(21) 



- "6 U " 4 iV^ " ivTiV^j + 8 " " iV:^^ + 12A^ iV:^^- ^^^^ 

Thus, we have obtained the QCD anomalous dimension and the ratio 
of average multiplicities in gluon and quark jets as functions of the QCD 
coupling constant 70 and of the number of active flavours in higher-order 
perturbative QCD. The higher order corrections in r are negative and its 
value becomes smaller. The above dependences are very mild that agrees 
with conclusions of [Q]. Besides, here we can insert the running coupling 70 
in (|13D and show dependence of 7 on Q directly ( see Fig.l ). To do that, 
we use ( as in 1^ ) the parameter Qo = 0.65A^ with = 175MeV for 
n/ = 5 in proportion 63 : 100 : 130 for n/ = 5 : 4 : 3, correspondingly. 
Integration of these values of 7 gives rise according to (|12|) to the following 
energy dependence of the average multiplicity 

/ \ ^ I a2 
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or in terms of 70 one obtains 



{nc) = 2/2^ exp[27o(?/)(l + a2^l{y))y] . (24) 

First term in the exponent corresponds to DLA-approximation with run- 
ning couphng constant. There is an important difference between the fixed 
and running couphng constant impact on the energy behaviour of average 
multiphcities. Even though the fixed couphng case provides faster asymp- 
totical increase of average multiphcity, in the energy range of interest up to 
y ~ 3?/^o the mean multiphcity for fixed coupling increases slower than for 
running coupling. ( Here yzo = 6.67 is the scale at which the value of the 
fixed coupling has been calculated). It is shown in Fig. 2. 

The physical origin of the effect is easy to understand because rather small 
value of the fixed coupling has been used throughout the whole evolution of 
the cascade while in real situation it should become larger and larger at 
later stages what produces more particles in case of running coupling at 
lower energies. However, the multiplicity is enlarged if more flavours become 
active. It compensates somehow the lower values of the coupling constant at 
high virtualities what should be taken into account when comparison with 
experiment is done. 



The preexponential term in expressions (23), (|2^) is also well known ( see, 
e.g., 1^, |ry] ) and stems from MLLA-approximation Q The next-to-next-to 
leading correction to the exponent calculated by us diminishes it and the 
energy increase of the average multiplicity becomes somewhat slower ( see 
Fig.3 ). 

In Table |l| we show the values of the ratio r for various n/ at Q = 91 GeV 
( mass ). They are very close to those obtained for the fixed coupling 
QCD and are about 20 percents smaller than the prediction of the lowest 
order approximation. We would like to point out that the 7Q-correction 



to r calculated above is several times larger than its value found in [10 



where the complicated diagram calculations have been done. The result 



of 1^ is reproduced by the first term in (0). Other terms correspond to 
the contribution of second derivatives in (|TUp,(|TTp. For convenience, the 
numerical values of shown in the Table also. 

Let us remind that one should combine these results with models for 
the conversion of partons into hadrons ( hadronization ) to get an access 
3 



Calculations in the framework of this approximation see, for example, in refs.|llL 12 
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to experimental data. That is why, in ref. 0, the results of Monte Carlo 
calculation according to Herwig model ||^ have been used. They show that 
at Z°-mass the ratio r at parton level exceeds its value at hadron level by 
about 40 percents. Being taken into account, this fact leads to agreement 
between theory and experiment |@] . Since our results are very similar to those 
of we do not foresee any contradiction between them and experimental 
data. 

To conclude, we have obtained the QCD anomalous dimension 7 and the 
ratio of average multiplicities in gluon and quark jets in higher-order pertur- 
bative QCD with running coupling constant. The value of the ratio agrees 
rather well with the one previously obtained in fixed coupling QCD by 
exact solution of the equations. Thus, it shows that the running property of 
the coupling constant is not as important for the ratio as proper treatment 
of higher-order terms including conservation laws but it is more noticeable 
in the energy dependence of mean multiplicity. Once again, it demonstrates 
that the values of r are much smaller than their lowest order counterparts 
and agree with experiment if the hadronization is properly taken into ac- 
count. Nevertheless, one may still worry that the hadronization stage gives 
rise to the most important factor reducing parton ratios to their hadron val- 
ues accessible in experiment. However, to get it in QCD one must solve the 
confinement problem that is out of our scope nowadays. 
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Figure Captions 



Fig. 1. The dependence of the anomalous dimension 7 on Q. Sohd hues are 
for 7 and dashed hues are for 70 . The numbers of active flavours are 
indicated near lines in all figures. 

Fig. 2. The y-dependence of the average multiphcity. Sohd hues are for 
running coupling and dashed lines are for fixed coupling. The arrow 
marks the y^o location. 

Fig. 3. The ^/-dependence of the average multiplicity for the running cou- 
pling. Solid lines arc for next-to-next-to leading approximation consid- 
ered in the present paper and dashed lines are for MLLA. 
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Table 1: 





r 


n 


r2 


7o 


ai 


02 


3 


1.84 


0.185 


0.426 


0.473 


0.280 


- 0.379 


4 


1.80 


0.191 


0.468 


0.481 


0.297 


- 0.339 


5 


1.77 


0.198 


0.510 


0.484 


0.314 


- 0.301 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/9406002vl 



This figure "figl-2.png" is available in "png" format from: 
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